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Abstract 



o 

■ The focus of the article is on analysis of skew-symmetric weight matrix functions for 
<^ . interfacial cracks in two dimensional anisotropic solids. It is shown that the Stroh formal- 

ism proves to be an efficient approach to this challenging task. Conventionally, the weight 
functions, both symmetric and skew-symmetric, can be identified as a non-trivial singu- 

■ lar solutions of the homogeneous boundary value problem for a solid with a crack. For a 
semi-infinite crack, the problem can be reduced to solving a matrix Wiener-Hopf functional 
equation. Instead, the Stroh matrix representation of displacements and tractions, combined 

, with a Riemann-Hilbert formulation, is used to obtain an algebraic eigenvalue problem, that 

is solved in a closed form. The proposed general method is applied to the case of a quasi- 
(— I , static semi-infinite crack propagation between two dissimilar orthotropic media: explicit 

expressions for the weight matrix functions are evaluated and then used in the computation 
of complex stress intensity factor corresponding to an asymmetric load acting on the crack 
faces. 



Keywords: Interfacial crack, Riemann-Hilbert problem, Stroh formalism. Weight func- 

^ , tions. Stress intensity factor. 

00 

! 1 Introduction 

■ Evaluation of coefficients in asymptotic representations of displacements and stress fields repre- 

ss! ! sents an important i s sue for addressing vector prob lems of crack propagation in elastic materials 



(jBercial-Velez et al.l . l2005l : iMishuris &: Kuhnl . |200l|) . The explicit derivation of weight functions 



is fundamental for the evaluation of stress intensity factors corresponding to a general distribu- 
^ . tion of forces acting on the crack faces, as well as for the calculation of higher order coefficients in 

\ the asymptotic expressions of the fields. The latter are used in asymptotic models of incremental 

^ ■ crack growth, and hence are essent ial for eva l uation of th e crack path and analysis of fracture 

stability. In the work by Bueckner ( Bueckner . 1985 . 19891 ) . weight functions for several types of 



cracks in homogeneous elastic media, both in two and three dimensions, have been defined. 

In this paper, the term "symmetric" load is associated to forces of the same magnitude ap- 
plied on both crack faces in opposite directions, while the load generated by forces acting on both 
crack faces in the same direction is called "skew-symmetric" or "anti-symmetric". For cracks 
in homogeneous elastic materials, in the two-dimensional setting the skew-symmetric loading 
does not contribute to stress intensity fa ctors, whereas it becomes relevant ari d it must to be 
accounted for in three-dimensional solids (Bueckner, IDSSl : Meade &: Keer . 19841 ). The situation 



is different when the crack is placed at the interface between two dissimilar elastic materials: 
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even for two-dimensional problems the skew- s ymnietric loads generate a non - zero contribution 



to str ess intensity factors (|Bercial-Velez et alj . l2005l : iLazarus &: Leblondl . Il998l : iPiccolroaz et al 



200^). In particular, for Mode III interfacial cracks the stress components do not oscillate, but 



a non-vanishing s kew-symme tric component of the weight function still has to be accounted 
(jPiccolroaz et alJ . hOQ^ . hoid ) . In the case of isotropic media, weight functions for semi-infinite 
cracks can be defined as singular non-trivial solutions of the homog eneous boundary problem 
with zero tractions on the crack faces but unbounded elastic energy (j Willis fc Movchanl . [l995l ). 
For interfacial cracks between dissimilar isotrop i c elastic media, the weight f unctions are well 
discussed in literature ( Lazarus Leblond . 1998 : Piccolroaz et al. . 20091 . 20ld ). The problem is 
generally reduced to a functioi ial equat i on of Wiener-Hopf type, and its solution gives the sym- 
metric weight function matrix ( Antipov . 19991 ). while the skew-symmetric component is obtained 
by the construction of the corresponding full-fiel d singular solution of the elasticity boundary 
value problem discussed in IPiccohoaz etaD (|2009l ^. Considering interfacial cracks in anisotropic 
dissimilar elastic media, instead, the results on skew-symmetric weight functions are not readily 
available. 

In this article we illustrate a general procedure for the calculation of symmetric and skew- 
symmetric weight functions matrices for semi-infinite two-dimensional interfacial crack problems 
in anisotropic elastic bi-materials. It is shown that the challenging analysis of the matrix func- 
tional Wiener-Hopf equation can be replaced by solving a matrix eigenvalue problem deduced 
via an equiv alent formulation, based on Stroh representation of displacement and stress fields 
(|Strohl . ll962l ). By means of this new approach, general expressions for weight functions matrices, 
valid for plane interfacial cracks problems between any anisotropic media, are obtained. This 
general result is reported and discussed in details in Section Q. 

Section Q illustrates the proposed method for the case of a semi-infinite two-dimensional 
stationary c rack between two dissimilar orthotropic materials under plane stress deformation 
(ISuol . ll990bl ). corresponding to a general non-symmetric load. Stroh representations for displace- 
ments and stress corresponding to this problem are explicitly calculated and used for deriving 
symmetric and skew-symmetric weight functions matrices. In Section ([5]) both symmetric and 
skew-symmetric weight functions are utilized together with Betti formula in order to evaluate 
complex stress intensity factor for an interfacial crack in orthotropic bi-material subject to an 
asymmetric loading configuration. In the particular case of isotropic media, the obtained result 
is consistent with the stress intensit y factor derived for the non-symmetric distribution of forces 
obtained by Piccolroaz et al. ( 20091 ). 

Finally, in Appendix A, the evaluated skew-symmetric weight function is compared to those 
calculated by the ful l field singular s o lution of the plane elasticity problem, following the ap- 
proach illustrated in IPiccohoaz elaP (200^). The perfect agreement detected between the ex- 
pressions derived by means of two alternative formulations represents an important benchmark- 
ing for the correctness of the obtained results. 



2 Interfacial cracks: preliminary results 

Here we introduce the main notations and the mathematical framework of the model. Consider 
a quasi-static advance of a semi-infinite plane crack between two dissimilar anisotropic elastic 
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materials with asymmetric loading applied to the crack faces, the geometry of the system is 
shown in Fig.([T]). 

The loading is defined via tractions acting upon the crack faces. Considering a Cartesian 
coordinate system with the origin at the crack tip, (see Fig.([T])), traction components behind 
the crack front are then defined as follows: 



4.(xi,0^ 



where p^{xi) are given functions. 



p - (xi) for xi < j = 1,2, 



(1) 



Since the load is assumed to be self-balanced, its resultant force and moment vectors are 
equal to zero. Moreover, we assume that forces are applied outside a neighbourhood of the 
crack tip and vanish at the infinity. The body forces are assumed to be zero. Symmetric and 
skew-symmetric parts of the loading are given by the following expressions: 



(Pj) i^i) = \ + Pj ixii) > [Pj]{xi) = P'j'ixi) - Pj (xi), j 



1,2, 



The solutions in form of functions which vanish at the infinity and possess finite elastic energy 
are sought. Expressions for the stress field and displacements for a semi-infinite interfacial c rack 
in anisotrop ic bi-materials have been obtained by means of Stroh formalism (jStrohl . Il962l ) by 
Suol (1990b). This approach to the physical problem of the crack, that is used for the derivation 
of weight function matrix in the paper, is reported in Section ()2.ip . 

In Section ()2.2p the weight functions are defined as non-trivial singular solutions of the 
homogeneous elasticity problem for interfacial crack s with zero tractions on the crack faces but 
unbounded elastic energy, following the approach of Piccolroaz et al. ( 20091 ). 

Section ()2.3p reports the application of the Betti formula to physical fields and weight func- 
tio ns and the derivat i on of t he fundamental identity in th e Fou rier space, discussed in details 



m 



Willis fc MovchanI (| 19951 ): IPiccohoaz k MishurisI (|201lh and IPiccolroaz et al.l (|2007l ). This 



integral identity will be used in the paper for the evaluation of coefficients in the asymptotics 
of the stress field s near the crack tip by means of a procedure based on weight functions theory 
(jPiccohoaz et al.l . l2009l ). 



The asymptotic representations of physical diplacements and stress f ields near the crack tip , 
including the stress iii t ensity factor as the coefficient of the first term ( Piccolroaz Sz MishurisI . 
201 ll : IPiccolroaz eTall . I2OO9I ) , are introduced in Section (12.41). 



2.1 Stroh formalism in analysis of interfacial cracks 

Physical displacements and stress fields for an intefacial crack between two d ifferent anisotropic 
mater i als can be derived by means of the Lekhnitskii or Stroh appro aches (jLekhnitskiil . Il963l : 
Strohl . Il962l : ISu 3. Il990bl ). Following the article by Suo (H, Il990bl ) we introduce the stress 



vectors tj = {(7ij,a2j)'^ together with the displacement u = {ui,U2)'^- The constitutive relations 
for both the elastic media oc cupying respe ctively the upper and lower half-planes can be written 
using the Stroh formulation (Stroh, 19621 ): 



ti 

t2 



Qu 1 + Ru 2 
R^u 1 + Tu i 



(3) 
(4) 
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Figure 1 : Geometry of the model. 



where the matrices Q,R and T depend on the material constants. A semi-infinite static crack 
placed at the interface between the two materials is considered, as it is illustrated in Fig.([T]). Ac- 
cordingly, the derivative of the displacement u^i(xi,X2) and the traction T{xi,X2) = t2(2;i,X2) 
can be written as 

r(xi,X2) = Bg(z)+Bg(z) (5) 

and 

ui(xi,X2) = Ag(z) + Ag(z), (6) 
where A and B are constant matrices, g(z) is an analytic vector function with compon ents 



gi(x^ + iijX2), and /ij are complex numbers with positive imaginary parts. According to ISuo 



|l990E), if gj{zj) IS an analytic function of Zj = xi + fJ,jX2 in the upper half-plane (or in the 



lower half-plane) for one nj, where iij is a complex number with positive imaginary parts, it is 
analytic for any /ij. On the basis of this property, here and in the text that follows, we reduce 
the analysis to a single complex variable. Th e connection between the elements of A, B and 
Q,R, T is given by the following relations (see Ting ( 19961 ). pages 170, 171): 



(Qik + {Rik + Rki)^^j + Tikn))Akj = (7) 

Bij = {Rki + fJ-jTik)Akj (8) 

Thus, each column of A is a non-trivial solution of the eigensystem ([7|), while the eigenvalues 
fj,j are roots of the characteristic equation: 

\Qik + (Rzk + Rki)f^j + Tikfi]\ = (9) 

In turn, for each of the two phases we introduce the Hermitian matrix Y = zAB~^, which will 
be used in the text below. Further in the text, we shall use the superscripts and to denote 
the quantities related to the upper and lower half-planes, respectively. 

The traction ([SJ can be expressed in both upper and lower half-planes in terms non-homogeneous 
Riemann-Hilbert problem: 

h+{z)-h-{z)=r{xi,X2), (10) 
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A crack advancing along the negative semi-axis xi < is considered, the free traction 
condition is imposed for xi < 0, while the continuity of the tractions at the interface is ass umed - 
ahead the crack front . The following equations are satisfied at the interface (for X2 — )• 0^) ( Sud . 
1990bh : 



h+(xi)+H Hh-(xi) 
h+(a;i) +H"^Hh-(xi) 



7~(xi) for xi > 
for xi < 



(11) 
(12) 



The detailed analysis of this problem is included in Suol ( 1990b ): Suo et al. (1992), and it 
shows that the stress and displacement fields do not have oscillation near the crack tip for 
the case when the matrix H = Y^^^ + Y^^^ is real, otherwise the stress and displacement are 
characterised by the oscillatory behaviour near the crack tip. The branch cut for the function 
h(z) is assumed to be along the crack line xi < 0. Assuming that the stress vanish at the infinity, 
the solution of the homogeneous Riemann-Hi lbert proble r n (1121) is h - (z) = wza"*"*^, where w is 
solution of the following eigenvalues problem (|Sud . Il990bl : ISuo et alifliii ): 



Hw 



„27re 



Hw. 



(13) 



The traction T ahead of the crack has the form: 



r(xi) 



:Re 



Kxfw 



(14) 



where K = Kj + iKjj is the complex stress intensity factor including both mode I and mode 
II contributions to the traction, e is the bi-material parameter (a real non-dimensional number 
measuring an aspect of elastic dissimilarity of the two materials), and w is the eigenvector 
obtained from ()13p . 

For the case of plane strain load, it has been shown in Suo (jSuol . Il990bl : ISuo et all Il992l l 
that the displacement jump [u] across a semi-infinite crack running along the negative semi-axis 
xi < is given in the form 



[u](xi) 



'2(-:ri)V/' (H + H 



vr 



cosh Tre 



Re 



Kj-xiY^w 
l + 2ie 



(15) 



2.2 Weight functions definition 



Following the theory developed by Willis and Movchan ([Willis &: Movchanl . Il995l ) , we define a 
vector function U = (f/i, f72)"^ as the singular solution of the elasticity problem with zero traction 
on the faces where the crack is placed along the positive semi-axis xi > 0. T he traces of these 



functions on the plane containing the crack are known as the weight functions (jPiccolroaz et al 



20091 ). and notations [U] and (U) will be used in the paper to denote symmetric and skew- 
symmeric weight functions respectively: 



[U](xi) = U(xi,X2 =0+) -U(xi,X2 = 0-) 
(U)(xi) = ^(U(xi, X2 = 0+) + U(xi, X2 = 0-)) 



(16) 
(17) 
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The traction vector 5] = (£i , S o)'^ associated to the s i ngular solutions is continuous on 
the plane containing the crack fsee Ipiccohoaz et al.l H, &)) and vanishes for xi < 
(homogeneous buondary conditions are imposed). In practice, for singular solutions we impose 
free traction condition for xi > and traction continuity at the interface for xi < 0: 



h+(xi) + H ^Hh"(xi) = for xi > 
h+(xi) +H"^Hh"(xi) = S(xi) for xi < 



(18) 
(19) 



It is important to note that the domain of singular solutions and then of weight functions is 
different from the domain of the physical solution, defined in previous Section, where the crack 
is placed along the negative semi-axis. For the singular solution, the branch cut for h(xi) is 



chosen to be along the line xi > 0. Using the representation h{z) 
the eigenvalues problem: 



Hv = e"2''^Hv. 



V2;2+*^ we reduce (fT9l) to 



(20) 



The bi-m a.terial matri x H is Hermitian, and the eigenvector corresponding to the eigenvalue —e 
is V = w ( Suol . 1990bl ). The singular traction is sought in the form: 



(-^i) 



Re(Cw(-xi)^ 



(21) 



where C = Cj + iCjj is a complex constant representing both mode I and mode II contributions 
to the traction (j2ip . as it is for the complex stress intensity factor K with regard to the physical 
traction (see eq. ()14p ). The shear and normal opening modes for plane stress and plane strain 
problems are coupled. Two linearly independent vec t ors S and associated w eight functions can 
be identified similar to Piccolroaz &: Mishurid (|201lh : IPiccolroaz et al.l (|2009l ^. 



2.3 Fundamental Betti identity 

As discussed above, U is discontinuous along the positive semi-axis xi > 0, whereas u is discon- 
tinuous along the positive semi-axis xi < 0. S is zero for xi > 0, whereas T is zero for xi < 0. 
Asymmetric loading are applied to the crack faces (see Fig. ([1]) ) , thus the physical traction acting 
at the interface on the entire xi axis can be written as: 



cr(xi,X2 = 0+) = p+(xi) T(xi), <t(xi,X2 = )=p {xi) + T{xi] 



(22) 



Ac cording to the approach illustrated in IWillis fc MovchanI \l99^ ) and iPiccolroaz et al.l 
(|2007l ). m order to derive explicit formulas for calculating the coefficients for the asymptotic 
of the stress fields near to the crack tip, we apply the Betti formula to the physical fields and to 
weight functions. The following rlations, respectively in the upper and in the lower-midplane, 
are obtained: 

/ |u'^(xi -xi,0+)7?,cr(xi,0+) - 5]^(xi -xi,0+)7?.u(xi,0+)|dxi = (23) 

Jx2=0+ ^ 

/ |u^(xi -xi,0")7?.cr(xi,0~) -5]'^(x'i -xi,0")7?.u(xi,0")|dxi = (24) 

Jx2=0- ^ -* 
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where TZ- is the rotation matrix: 



-1 
1 



Subtracting the (j24p from the (j23p and using (j22p . we derive the following integral identity: 
/ {u'^ix'i -xi,0+)7^p+(3;i) +U^(x'i - a;i,0+)7^r(xi) - U^(a;'-L - xi, 0") - 7^p-(a;i)- 

JX2=0 ^ 

-U^(xi -xi,0")7?.T(xi) - -xi,0+)7?.u(xi,0+) + 5]^(x'i - xi, 0")7?.u(xi, 0")]} dxi = 

Referring to (fT6]) and (fT7|) we deduce 

{[U]^(x; - xi,0+)7?.r(xi) - S^(x; - xi,0+)7^[u](xi)} = 



L 



X2=0 



• / ^[U]^(x;-Xl,0+)7^.(p)(xl)-S^(x;-Xl,0+)7^[p](xl)|. 



(25) 



Let us define the Fourier transform of skew-symmetric weight function and associated trac- 
tion 5] respect to xi: 



/■+00 _ ^0 

[U]+(e)=/ [U](xi)e^«-Mxi, t (0= 5](xi)e*«-Mxi, 

JO J-oo 



(26) 



Where the superscript -|- indicates that [U] (^) is analytical in the upper half-plane (Im^ E 
(0, +oo)) and the superscript — indicates that XI (^) is analytical in the lower half-plane (Im^ G 
(—00,0)). The physical traction and the jump function are defined in such a way that the 
transforms 7" {£,) and [u]~(^) are analytic in the upper and in the lower half-plane, respectively. 
Applying the Fourier transform to t he (1251), by means of convolut i ons properties, the foll owing, 
relations, valid for ^ E R, is derived (jPiccolroaz &: Mishurid . 1201 ll : IPiccolroaz et al.l . l2007l ) : 



[U]+^7^(p) - (U)+'^ 7^[p], (en 



(27) 



This identity, obtained by Willis &: Movchan ( 19951 ) and Piccolroaz et al. ( 2007 ). relates trans- 
forms of the physical solutions to transforms of the weight functions, and it is used for the 
evaluation of the stress identity factors. 



2.4 Asymptotic representation of the fields and stress intensity factors 

The physical traction ()14p and the displacement jump across the crack face (|15p for xi — ?■ 0, can 
be written in the matrix form, as follows: 



r(xi) 



2V27r 



T(xi)K + 



X 



[u](xi) = ^-^^U{x{)^ + 



2V27r 

(-^i) 



T(xi)A + 



"1 



2v/27r 



T(xi)B + 0(xf) 



27r 



'-U{xi)A + 



-Z^(xi)B + 0((-xi)2; 



(28) 



(29) 
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Where K = (K, K),A = {A, ^4), B = (B, B), and A and B are higher order coefficients, defined 
in the same way of the stress intensity factor: A = Aj -\- iAjj,B = Bj + iBjj. Matrices T(xi) 
and U{xi) are given by: 



T(xi) = 2(^ 



-1 1 



2(H + H) (wi-xiY^ W{-xi] 



cosh Tie 



l + 2ie 



1 - 2ie 



Th e evaluation of coe f ficient s in the asymptotic formulae (j28p and (j29p is included in lPiccolroaz &: Mishuris 
(|2nilh : IPiccolroaz eFID (120091 ) :or interfacial cracks between dissimilar isotropic media, by 
means of a general integral formula involving symmetric and skew-symmetric weight functions 
matrices. Here this approach is extended to anisotropic bi-materials: in the next Section general 
expressions for Fourier trasforms of symmetric and skew-symmetric weight functions are derived 
using the Stroh formulation, as defined in Section (|2.ip . while in Section ([¥!) these ex pressions 
will be specialized to the case of a crack between two orthotropic materials (ISuo . 1990b) and the 



gener al formula obtained by Piccolroaz et al. (jPiccolroaz fc Mishurisl . 12011 
20091 ) will be applied for calculating stress intensity factor for this case. 



Piccolroaz et al. 



3 Symmetric and skew-symmetric weight functions for anisotropic 
bi-materials 

Here we derive general expressions for symmetric and skew-symmetric weight function defined 
in Section (|2.2p as the jump [U] and the average (U) of the singular solution for a semi-infinite 
interfacial crack with free-traction conditions. 

Let us introduce the vector function h(2:), solution of the Riemann-Hilbert problem for the 
physical traction ([TO]) , such that: 



h(^) 



Bg(z), Im z > 0, 
— Bg(z), Im z < 0. 



According to the Plemelj formula, h(2;) can be written as follows: 

r(r) 



1 /■'^ 



T — z 



-dr, 



where the integral is understood in the principal value sense. Taking the Fourier transform 
respect to xi, for X2 = 0^, at the interface, we obtain: 

/oo I 
-oo 

where H is the Heaviside function, and the traction T is defined in such away that its transform 
is analytic in the upper half-plane. To derive the above representation we used the following 
relations 

/oo pi£,{xi±iO) 
-oo T - [xi ± ^0) 
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Hence the Fourier transforms of g and g at the interface are obtained in the form 

m = B-ih+(o = i/(-OB-ir"^(o, i G R, 

and 

t(e) = -B"'h"(e) = H{i)B-'t^{0, i G R. 

By applying the Fourier transform to ^ for the derivative of the physical displacements, we 
deduce: 

-ieu(e, X2 = o±) = Am + At(e), e e r 

Hence, the Fourier transform of the displacements on the boundary of the upper half-plane is 
given by: 

u(e,x2 = 0+) = ^{i?(-0AB-i + F(OA B"'} r^(e) 

By expressing the Heaviside function in the form: 

/7(±e) = ^(i±sign i) 

we obtain: 




Following the same pattern of the derivation as above, we derive the Fourier transform of the 
physical displacement on the boundary of the lower half-plane: 




Next, we replace u in the above text by the singular solution U for a semi-infinite interfacial 
crack. Correspondingly, the vector of tractions T is replaced by XI, which is the traction vector 
corresponding to the singular solution U. The Fourier trasforms of singular displacements on 
the boundary can then be derived: 




According to definitions (fT6]) and (fT7|) . we take the difference and the average of (j32]) and (p3|) . 
and we derive that: 

[U]+(e) = ^^i sign(e) Im(Y(i) - Y^^)) - Re(Y(i) + Y^^^)}±- [0. (34) 
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and 



(u>(0 



^[i sign(e) Im(Y(i) + Y(2)) - Re{Y^^^ - Y^^^)}± (0, ^ € R- 



m 



(35) 



We note that (l3^ is the functio nal equation of the Wiener-Hopf type, similar to the one studied 



m 



Piccolroaz et al.l (|2009l . l2007l ) and lPiccolroaz k MishurisI (|201ll ) for the case of isotropic media. 
Using the continuity of tractions S across the interface, together with ()34p and ([35 



express the Fourier transform of the skew-symmetric weight function (U) in the form: 

(u)(e) = ^[u]+(e) + >5]"(0, eeR, 

where the diagonal matrix A. and the off-diagonal matrix B are 



we 



(36) 



A=^ Re(Y« - Y(2))(Re(Y« + y(2)); ' 



and 



B 



- Im(Y(i) + Y(2)) - A Im(Y(i) - Y(2)). 



The formulae (j34p and (j35p give expressions for the Fourier transform of the symmetric 
and skew-symmetric weight functions for interfacial cracks in anisotropic bi-materials in terms 
of the transformed singular traction S. These expressions are compared to those obtained by 
the construction of the full-field singular solution for the elasticity problem in an half-plane 
(jPiccolroaz et all bond l in Appendix A. As expected, the weight functions derived via two dif- 
ferent approaches agree. Making the inversion of the (j34p and the (j35p the explicit representation 
of the weight functions can be obtained and the complex stress intensity factors associated to 
an arbitrary system of forces car i be evaluated by mea ns of integral formulas derived by the 
application of the Betty identity ( Piccolroaz et al. . 20091 ). 



4 Interfacial crack in orthotropic bi-materials 

In this section explicit expressions for the Fourier transforms of the symmetric and skew- 
symmetric weight functions for an interfacial crack in orthotropic bi-materials are derived using 
equations and (f35|) . 



4.1 Stroh representation for orthotropic bi-materials 

For the case of orthotropic two-dimensional media, the matrices Q,R and T introduced in 
section ()2.ip are given by: 



Q 



Cll 



066/' 



R 



C66 
C22 



Ci2 

,C66 

Where Cij are the elements of the materials stifness matrix, which can be expressed in function 
of the elements of the compliance matrix Sij (Suo, 1990b): 



Cll 
C22 





S22 


C12 = 




S12 


*12 


- SIIS22 


*12 


- S11S22 




fill 


C66 = 


1 




■^12 


- SIIS22 


S66 





(37) 
(38) 
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The characteristic equation (jlOp becomes: 

Sii/i^ + (2si2 + S66)/i^ + S22 = 



(39) 



The matrices A and B, defined by relations (l8|) and ([8]), are equivalent to those provided by 
alternative Lekhnitskii formulation (jLekhnitsk ii. 1963), more precisely the Lekhnitskii approach 
gives a specially normalized eigenvector matrix A, and the characteris t ics eq uation derived using 
the Lekhnitskii formalism is identical to the ()39p ( Suol . 1990b : Ting . 19961 ). The relation ships 
betwe en the two alternative formulations are derived and reported in details in reference iHwu 
(|l993l ). where the elements of A and B in the Stroh representation are reported in function of 
the coefficient proposed by Lekhnits kii. H e re we write the Stroh matrices A and B using this 
particular normalization, reported in iHwu d 19931 ): 



y^^(s22-SllMi) ^J'^^i^22^-Siu4) 



\ 



S12M2 + - 



V V7^(''22-SllM^) 
i£3_ 



(40) 



J^(s22-SllMt) / 
-^2 



B 



,J^^{S^22^^I^U4) y^7^(«22-SllM|) 
1 1 



(41) 

V \/w(^22-Siim4) y'-l-(s22-Sii/x|) j 

Where //i and ^2 are the two roots of the characteristics equation with positive imaginary part. 
The Hermitian matrix Y evaluated using (|4Up and (|4ip is: 



iAB 



siilm(/zi + ;U2) -i(sii/Ui/U2 - S12) 



y ^^Sll/Xl/i2 - Sl2j -S22lni(^ 

Following the notation introduced by Suo, (|Suol . Il990bl : ICuota et al.l . Il992l ). we define two 
adimentional parameters measuring the material anisotropy: 




(42) 



S22 



_ 1 2£i2_+£66 
2 V^llS22 



If A = 1 the material has trasversely cubic symmetry, while if A = p = 1 the material is 
trasversely isotropic. The positive definetess of the strain energy density requires that: 

A > 0, -1 < p < +00, 

The characteristics equation ()39p then becomes: 

A/ + 2pA5/i2 + 1 = 

The roots with positive imaginary parts are: 

jii = iX^*{n + m), /i2 = iA~4 (^ — rn), for 1 < p < +00, 
pi = X^i {in + m), fi2 = X^^{in — m), for — 1 < p < 1, 



(43) 



P2 = «A 4, 



for 
for 



1, 
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Where: 



n 



Using this notation the matrix Y becomes: 

\-i((siiS22)2 + S12) 2nA"4(snS22)2 / 
The bi- material matrix H for two orthotropic materials ( Suo . 1990bl : Gupta et al. . 19921 ) is: 



H 



-if3^HuH22 
if3\/HiiH22 H22 



(45) 



Where: 



Hii = [2n\l{siiS22)^^^ + [2n\^{siiS22)H^\ 
H22 = [2nA-3(siiS22)5](^) + [2nA~^(snS22)^]^^\ 

NH11H22 = [((S11S22)^ + Si2)](') - [((S11S22)^ + Si2)]^^\ 

P is the generalization of one of the Dundurs parameters ( Dundurs . 19691 ). and is connected to 
the bi-material oscillatory index e by the relation: 

The eigenvector w of the eigensystem (jl3p . associated with the component of displacement jump 
(1151) and of the traction ahead of the crack (1141). is assumed in the normalized form: 



1 1 H 



T 



W 



11 



2 ' 2 \ H' 



22 



(46) 



4.2 Weight functions for in-plane deformations 

Relations ()34p and ()35p gives general expressions for the Fourier transform of the symmetric and 
skew-symmetric weight function in function of the transform of the traction SI. The singular 
traction on the boundary 5] for an interfacial crack between two generic anisotropic materials 
is given by ()2ip . Here we report this expression: 



^iL^Re fC7w(-xi: 



(47) 



In order to study the propagation of a crack between two orthotropic media we assume for w 
the expression (fi6|) . As just anticipated, mode I and mode II are coupled, and associated to one 
single complex constant C . Therefore the spaces of singular tractions and of singular solutions U 
are two-dimensio nal linear spaces, and two linearly independen t vecto rs 5] and weight functions 
must be defined ( Piccolroaz et al. . 20091 : Piccolroaz &: Mishurid . 2011 ). Assuming = 1, C// = 
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0, and Cj = 0, Cu = 1, and substituing explicit components of w, the two independent traction 
vectors are defined: 



i-xi)- 



2V27r 



2V27r 



>h/fT((-^ir-(-^i)""). 



(48) 



(49) 



Applying the Fourier trasform to (fl8|) and (fM]) . the result is: 



e 



l + 4e2 



7^ 



( 



l + 4e2 

Where Im^_ E (— oo,0), eo = e^^ and: 



V V ^^22 



(-i + ie 



(50) 



(51) 



± (l + «)\/vr 
eo = 2 , - ' 



2r 



± ie 



Substituing the (j50|) and (j50|) into (j34|) and (j35|) , and expressing the matrices in term of elements 
of H, we derive the Fourier transforms of the two independent symmetric and skew symmetrix 
weight functions: 



21+ ■ 



VH11H22 



( 



2J [^2J 



#^ i/3sign(e) 





V-i/3sign(0 v^f^ 



i7sign(^) ^2\[^ 









Pi]-/ 






\[^\]- 


Pi]-. 



(52) 



(53) 



Where the following Dundurs-like parameters have been defined: 

[2nA3(siiS22)^]^^^ - [2nAi(siiS22)^](^) 



7 



H 



11 



[2nA-3(siiS22)^]^^^ - [2nA-i(siiS22)^](^) 



22 



[((S11S22)^ + gl2)]^^^ + [((^11^22)^ + 512)]^^^ 
y/HiiH22 
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The superscript + indicates that U^(^) is analytic in the upper half-plane (Im^ G (0, +00)). 
The expression (j52p for the transform of the symmetric weight functions, as just anticipated 
in previous session, is a Wie ner-Hopf type equation simi l a,r to the one solved for inte rfacial 
crack in isotropic materials in iPiccolroaz &: Mishurid (l201lh : IPiccolroaz et alJ dioOTl .boO^'). The 
skew-symmetric part of the weight function can also be decomposed by means of relation ([53]): 

((]) = A[IJ]+ + jBt' , (54) 

Where the matrices A and 15 for orthotropic bi-materials using our notation are given by: 



B 



1 + P5i 

+ 



The explicit expressions for the weight functions are obtained by the inversion of the derived 
transforms. The symmetric weight function [U](xi) is equal to zero for xi < 0, while for xi > 
it is given by: 



Ul\ (xi) 
(^1) 



Hux^ 



2c+c-^/2^(l -F4e2) 
_ i_ 



(/3-l)(-- + ie)x-^ + (/3 + l) 



1 



2c+c-\/2^(l -F4e2) 



2c+c-^/2^(l + 4e2) 
VThlH^x^^' 



(/?-] 



+ ie ] x 



2c+c-V27r(l -F4e2) 



(/3 + 1 



(/3-l)(-- + ie)x-^ + (/3 + l) 



le I xf 



+ le ] X 



l£ X 



---ie]xf 



(55) 



The skew-symmetric weight function (U)(xi), is equal to ^[U](x) for xi > 0, while for xi < 
it is given by: 



{Ul){x,) 



iffll(7 + /j^l)(-Xi)-: 

4c+c-^/2^(l + 4e2) 



1 



+ ie 



+ 



VHnH22{7 + PS2){-xi] 
4c+c-\/2^(l -F4e2) 

Hi^{j + P6i){-xi)-^ 



+ ie 



+ ie 



-XlJ 



4c+c-\/2^(l + 4e2) 

WHuH22h + /362){-xi)-"2 ( (I 
4c+c-^/2^(l + 4e2) ^ 1 " 



(-^1)- 



1 

2 " 
- + 



.e)eg(-xi)-| 
ie^ eo(-xi)* 
ie]el{-x,r] 



-2-i£]eo{-xiY 



(56) 
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4.3 Stress-intensity factor for orthotropic bi-materials 

The symmetric and skew-symmetric weight functions are now used for evaluating the complex 
stress facto r for a loading correspondent t o an arbitrary system o f force s, following the procedure 
outlined in IPiccohoaz etaD (|2007l . |2009|) and IPiccolroaz et alJ (|20ld l. In Section [M we have 
shown that complex stress identity factor represents the coefficient of the first order asymptotic 
term of physical traction field bnear the crack tip. Remembering equation (|28p . for X2 — ?• 0"*" 
the traction becomes: 



r(xi) 



2V27r 



T(xi)K + 



X 



2V27r 



T(xi)A + 



3 

_xl_ 
2x/2¥ 



T{xi)B + 0{xl 



(57) 



For orthotropic bi-materials, considering the components of w given by (|46p . the matrix T(xi) 
is: 



-IX 



IXi 
^22-^1 



T(xi) 

The Fourier transform of the ([FT]) , as ^ — )• oo and Im^ G (0, -|-oo), is: 



t 2 2 £!r 2 

^+ t(e+)K + %^t(e+)A + %^t(e+)B + o(r 



Where: 



4^2 



(58) 



t(e4 



H22 c+eo 



H22 c~ 



7 



Using this expression, the explicit transforms of the traction 5] and of the symmetric 
weight functions matrix [U]^ derived in previous S ection, and evalu ating [u]~ , the asymptotic 
expansions for the members of Betti identity ()27p (jPiccolroaz et al ], I2OO9) are derived: 



-T^^+ _ ^-i_A4iK + C'^M2A + O(r^) for Im^ € (0, +00) 
1: ^n[u]- = C^MiK + C'^M2A + 0{C^) for Im^G(-oo,0) 

The explicit form for the matrix A4i is: 

^_(£^i)(i^ e2(/3 + l)(l + 2ze) 



(59) 
(60) 



Ml 



Hn 



4c+c-(l + 4e4) 



i{l3-l){l-2i£) 
? 



ieg(/3 + l)(l + 2ie) 



(61) 



Now, following Piccolroaz et al. ( 20091 ). we rewrite the Fourier transform of the Betti identity 
271 ) in terms of a Riemann-Hilbert problem: 



*+(e) - *-(o = -[u]+^(e)7^(p)(o - {vr'xcmmo, e e r 



(62) 
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Where ^(C)) according to the Plemelj formula, is: 

*(t) 



Then the solution of the (1621) is 



1 /-oo 



[t]+^7lf = Im^ G (0, +oo); 

S"^72.[u]" = Im^G(-oo,0); 

From these expressions asymptotic estimates can be extracted. For ^ — )• oo we can expand 
Plemelj 's formula, considering only the first term we have: 

= Si 737^- = + (64) 

Substituing this expression and the expansions ([Mj) and in comparing the correspond- 
ing terms, and considering only the first order, the following general formula for the complex 
stress intensity factor is obtained: 



K 



1 r 1 

—Air' / {[U]+^(r)7?:(p)(r) + (U)+^(r)7?:[p](r)}rfr (65) 



Since we have explicit expressions for matrix Ali and for the Fourier transforms of both sym- 
metric and skew-symmetric weight functions, now using the formula (j65|) we can evaluate the 
complex stress intensity factor associated to an arbitrary loading for interfacial cracks in or- 
thotropic bi-materials. In the next section an illustrative example of computation of K by 
means of the (|65]l is reported. 



5 An illustrative example 

In this section we present an illustrative example of computation of complex stress inensity 
factor K for an interfacial crack loaded by a simple asymmetric force system in orthotropic bi- 
materials. The considered force system is illustrated in Fig. ([2]): the loading consists in a point 
force F acting upon the upper crack face at a distance a behind the crack tip and two point forces 
F 12 acting upon the lower crack face at a distance a — 6/2 and a -|- 5/2, respec tively, behind the 



crac k tip. The loading can be expressed in terms of the Dirac delta function (jPiccolroaz et al 
[20091): 

_ F F 

p+(xi) = -F(5(xi -I- a), p~(xi) = -— ^(xi -I- a -I- 6) - — 5(xi -I- a - 6), (66) 

As it is shown in Fig. ([2]), the loading can be decomposed into symmetric and skew-symmetric 
part: 

F F F 

{p){xi) = -—6{xi + a) - -^S{xi +a + b) - -^^{xi + a - b) 

F F 
\p] (xi) = -F5{xi + a) + -^S{xi +a + b) + -^^{xi + a - b) (67) 
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F/2 



Asymmetric loading 



Symmetric part 



Anti-symmetric part 



Figure 2: "Three point" loading of an interfacial crack. 



The Fourier transform of the symmetric and skew symmetric part of the loading are given by: 

mo 

[p] iO = -Fe-'^" + |^e-*«('^+^) + |^e-*«('^-'') (68) 



Using these expressions and the exphcit transforms of the symmetric and skew-symmetric weight 
functions in orthotropic media, ()52p and (j54p . both symmetric and anti-symmetric part of the 
complex stress intensity factor K = + corresponding to this loading system have been 
evaluated by means of the integral formula (j65p : 



22 /2 __i 



1-/3' y i^ii V vr 
(1-/3)2 



a 2 



i + i(l-6/a)- 



+ ^(l + 6/a)-5-^- 



In order to study the behaviour of these symmetric and anti-symmetric contributions to the 
stress in tensity facto rs in function of b/a, the following non-dimensional parameters have been 
defined (|s^ . 1 1990^ 1: 



[(S11S22)^](') 



Sl2 



(1) 



Sl2 



(S11S22) 



[(S11S22)2](^) L(S11S22 

Then we can exprime Hu, i^22) <^2 and 7 in function of these parameters: 

H22 _ [2nA-i](^) + [2nA-i](2)^ 
Hu ~ [2nA3](i) + [2nAi](2)$ 



(2) 
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^ _ [2n\--4^) - [2nA-i](^)$ 

^ ~ [2nA-3](i) + [2nA-i](2)$ 

^ [l + 9]W-[l + 9](^)^ 
^ ^$[1 + G](i) - [i + e](2) 

It is important to note that the Dundurs parameter 7, associated to the skew-symmetric part 
of the loading, depends o n (3. Moreov er, oscillations of the stress and displacement fields are 
excluded for /3 = 0, e = (jSuol . |l99o3). In this case also 7 = 0. 

The complex stress intensity factor has been computed for /j^^) = 0.74, l/A^^) = 1,9(1) = 
1/2, p(2) = 4.91,$ = 6.4 (here we have considered that material (1) is alluminium and material 
(2) is boron (ISuol . Il990al )). G^^^ = 1/2,9(2) = 2 and five different values of the Dundurs oscil- 
lation parameter (5 = {—1/4,-1/2,0,1/4,1/2}. The values have been normalized multiplying 
by a'^F~\ and plotted in Fig.([3|) in function of the ratio b/a. The symmetric stress intensity 
factor is reported on the left of the figure, while the skew-symmetric is on the right, the real 
part is reported on the top while the imaginary is on the bottom. Observing the figure, we note 
that both the real and the imaginary part of the skew-symmentric stress intensity factor are 
zero for b/a = 0, as expected, since for 6 = the loading is symmetric. As we increase b/a, 
the skew-symmetric contribution to the loading become more relevant, and the skew symmetric 
stress intensity factor correspondingly increases. As we can see, in the case without oscillation, 
corresponding to /? = 0, both the imaginary parts Kjj and Kf'i vanish and the stress intensity 
factor becomes real. Both the symmetric and the skew-symmetric stress intensity factors diverge 
for 6/a — )• 0, because a point force is approaching the crack tip. 

In order to characterize the magnitude of the skew-symmetric stress intensity factor respect 
to the symmetric stress intensity factor, the ratio Kf /Kf is evaluated as a function of b/a in 
Fig-Q- We observe that as h/a increases, Kf may reach the 40% of , as a consequence 
we can say that the contribution of the skew-symmetric part of the loading is not negligible, 
and needs to be taken into account in perturbative analysis of interfacial cracks between two 
dissi milar anisotropic e l astic r nateri als subject to asymmetric forces systems applied on the crack 
faces (|Piccolroaz et al.1 . l2O10l . boOfll V 
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Figure 3: Symmetric and anti-symmetric stress intensity factors as functions of h/a computed for p'^^ = 
0.74, l/A'i^ = l,e'^^ = 1/2, p(2) = 4.91,1/A(2' = 14.3, 9^^' = 2,$ = 6.4 and different values of the Dundurs 
parameter P: P = -1/2, P = -1/4, /3 = 0, ^ = 1/4, /3 = 1/2. 
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0.2 0.4 0.6 0.8 1 

bla 

Figure 4: Ratio Kt/Kf computed for p'^' = 0.74, l/A^i' = = 1/2, p^^^ = 4.91, l/A*^' = 14.3, 8'^' = 

2, $ = 6.4 and different values of the Dundurs parameter ^: /3 = -1/2, 13 = -1/4, ^ = 0, /3 = 1/4, /3 = 1/2 
reported in function of h/a. 



6 Conclusions 



The developed general approach for the derivation of symmetric and skew-symmetric weight 
function matrices for interfacial plane cracks between dissimilar anisotropic materials, based on 
Stroh formulation of displacements and stress fields, have been discussed in details and tested by 
means of the application to the case of a crack placed at the interface between two orthotropic 
materials under plane stress. The skew-symmetric weight functions obtained by means of the 
proposed method have been compared to those obtained for the same proble m by the construc- 



tion of the full-field singular solution of the elasticity problem in a half-plane (jPiccolroaz et al 
200^), the comparison between the two different solutions is reported in Appendix A. The per- 



fect equality detected between the expressions derived by means of two distinct approaches is 
an important ulterior proof for the obtained results. 

Since the proposed Stroh repres entatiq i i is valid for stat i onary and steady-state elasticity 
problems in many anisotropic media (IStrohl . Il962l : ISud . ll99nbl : lTing] . ll99fil '). it can be utilized for 
evaluating explicit weight functions for plane interfacial cr acks in several kind of m aterials (Stroh 
analysis have bee n proposed for example in quasi- crvstals (iRadi fc Marianol . l20ld ^ , piezoelectr 



(jSuo et al.l . I1992I ) and poroelastics media (jGautier et al.l . I2OIII ) . The derived weight functions 



tries 



can be used in many i niportant applications: in perturbative expansions for growing cracks or 



wavy cracks problems ( Piccolroaz et al. . 20071 ). and in the computation of the stress intensity 



factor for non-symmetric self-balanced load generated by a system of point-forces applied on 
the crack faces. An example of stress intensity factor evaluation for an asymmetric loading 
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is reported in Section ([5]): weight functions matrices obtained for orthotropic bi- materials have 
been used in the computation, and the results show that the contribution of the skew-symmetric 
part of the loading is not negligible and must be considered in the asymptotic expressions of the 
stress near the crack tip, as it has already been demonstrated for the case of isotropic media 
(jPiccolroaz et alJ . bood boid ^ . 
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Appendix A 

In this appendix, the Fourier transforms of the singular solution of the interfacial crack problem 
between two dissimilar orthotropic materials are derived by solving a boundary value problem 
for a semi-infinite half-plane subjected to tract i on bo undary conditions at its boundary, following 
the procedure illustrated in iPiccolroaz et al. I (|2009l l. The derived expressions for the singular 



displacements and for the symmetric and skew-symmetric weight functions are compared to those 
obtained in Section ([3|) by means of the direct solution of the Riemann-Hilbert problem (jlOp . 
The perfect agreement detected between the expressions derived using two different approaches 
is an important test for the obtained results. 

Initially, we consider the lower half-plane, denoted in the article by the superscript 
Introducing the Fourier transform of the stresses respect to the variable xi, we consider the 
component (T22 as the primary unknown function, so that the plane strain elasticity problem for 
the orthotropic material ^"^^ is reduce to the following ordinary differential equation: 



^(2) - _"" _ c2 f (2) , 9 (2)\ . -" M (2) . _ _ 
'11 ''^22 ? [^66 ' ■^^12 J '^22 ' ^ '^22 "^22 ~ ^ 



Where a prime denotes the derivatives respect to X2. The characteristic equation associated to 
(IMI) is: 

[-^^¥4? - e (4? + 2^2^ [-^'^f + e^4? = 0, e G R 

Introducing r/^^^ = uj^"^^ /|^| and using the same notation of section^ this characteristic equation 
bcconiGS' 

A(2)[^(2)]4 _ 2p(2)(A(2))|[^(2)]2 + 1^0 (70) 

Assuming that this equation possesses four distinct roots, (p*^^^ 7^ 1), the general solution of the 
plane strain elasticity problem in the lower half-plane is: 
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The fourier transform on the displacements components are: 



u 



1 ~ t I ^11 '^11 + '^12 ^22 j ' ""2 — ^2 I '^11 "^11 ^12 ^22 ) ' 



(72) 



Where only the two eigenvalues with positive real part (Re^^^jy > 0), shuch that the stresses 
vanish at the infinity {a~j — )• for X2 — )• — oo), has been accounted. Remembering the 
conditions for having positive definetess of the strain energy density introduced in section [H the 
two eigenvalues with positve real part become: 



m 



(2) 



(2) 



A 4 (ri + m) 



A 4 (n + im) 



(2) 



(2) 



(2) (2) 



A 4 (n — m) 
A~4 (n — im) 



(2) 



for 1 < p^"^^ < +00 



(2) 



for 



1 < < 1 



From this form it is straightforward to note that these eigenvalues can be expressed in function 
of Stroh eigenvalues introduced in second section by means of the relation: 



Vi 



(2) 



■ (2) 



(2) 
V2 



■ (2) 
-^^2 ' 



(73) 



In order to derive the weight functions, we need to evaluate explicit expressions for the singular 
displacements utilizing the (|72|) . The boundary conditions along the boundary X2 = 0~ are 
defined as follows: 

<5"22(?'2;2 = 0") = E22(0> 0-2i{^,X2 = 0~) = t^^{^), ^GR 



where S22 components of the singular traction defined in section UJ (equations 

and (IMl)). It follows that: 



<5-22(C,a;2 = 0" 



-^slg: 



:n(e)K),f^ + 4%r)=^2i(0, 



and thus: 



A 



(2) 



^?^^22 - «sign(^)S2i 



(2) (2) 

V2 -Vi 



(2) _ isign(0S2i -i]f^T.^2 

(2) (2) 



A 



The Fourier transforms of the singular displacements fields are then: 



(74) 
(75) 



Ml (C,X2) 
1 



Vi 



sign(^) ( sjj - sS? 



Vi 



(2) 



^21 +W2 



(2) L(2) 
'12 



,(2) 
'11 



Vi 



(2) 



^22 



+ 



sign(0 ( sf} 



V2 



(2) 



■S12 ) ^21 + 'iVl' ( S 



(2) ( „(2) 
^11 



V2 



(2) 



*12 



^22 



\c\ (2) 



(76) 
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U2 {C,X2) 



iU?-r,f^) 



I s 



,(2) 
'11 



(2) 



1 2 



(2)_ (2)\ (2) 
*12 *66 ) '/I ^21^^ 



+sign(?) ( 45 + 45 - 4? 



(2) 



(2)J2)^. 
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(2) , J2) (2) 



12 ^ ''66 ''11 



■ (2)1 2 



+sign(0 ( 



(2) 



^(2) (2)\ (2) (2) 



^12 ^66 ^1 % 



^22 



(77) 



For the upper half-plane, we find the same expressio ns, subject to replacing |^| with — 
and the superscript (2) with (1) Piccolroaz et al. ( 20091 ). From equations ([75]) and ([77]) and 
their corresponding expressions on the upper half-plane, we can derive the traces of the singular 
displacements transforms on the plane containing the crack: 

/ [(??2 + ??l)sil]^^^ .[Si2 + Siir]i7]2]^^^\ 

{ 1 j"^ 

/ . [(sii + See) - {m + m + ??i??2)sii]^^^ [siir?ir/2(r/2 + m)]^^^ \ ^- 









{^,X2 


= 0+ 


lit 


(6 


= 




= 0+ 




(0 


= Ui 


{^,X2 


= 0" 


U2 


(0 


= U2 


{^,X2 


= 0" 



iim + m)'5ll]^^^ . [Sl2 + Slim^2]^^^ 

— ,^ 



1^1 



(sii + See) - (??2 + m + ??l??2)sil]^^^ [sil??l?/2(??2 + ??l)]^^^ 



V ^ 

Using the relation ([73|) between 7/1^2 and the Stroh eigenvalues: 

rji = -im, r]2 = -iH2, 
and considering the following relations 

2,2 / 2 I 2 

Vl+V2 = -(Ml +^2 



2,2 / 2 , 2\ 2si2 + See 22 22 -522 

' - - ■■ ^- , mv2 = ^^l^J■2 = — ; 

Sll Sll 



we deduce the expressions for singular displacements along the axes of propagation of the crack 
(x2 = 0), as follows: 







= ut 


(.^,X2 


= 0+ 


u'2 


(0 


= ^2 


{^,X2 


= 0+ 




(0 


= 


{^,X2 


= 0" 


U2 


iO 


= U2 


{^,X2 


= 0" 



[Im(/Ui + /i2)sii](^) . [Sii^l/i2 - S12 



(1)' 



. [siUn/^2 - Sl2]^^^ J_ 

■[Im(/ii +/i2)sii](2) 











Im ( — 




VAii 





[Sii/ii^2 - S12] 



(2)^ 



[Sll^l/J2 - S12 



-Sl2l(2) _J_ 







(2)\ 


Im ( — 






V/"i 


1^2/ . 





(78) 
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These expressions can be written in the same form of the physical displacements (|30p and ()3ip : 
U+(e) = {^(Y« - Y«) - ^(Y« + y(^))}e (79) 
U-(e) = {^(Y(^) - y(^^) + ^(Y(^) + Y(^))}S (80) 



where the hermitian matrices Y^^^ and Y^^-* possess exactly the same form ()42p . evaluat ed by 
specializing the Stroh formalism to the case of a two-dimensional orthotropic material (jSud . 



1990bD: 



Y(1)'(2) 



[snIm(^i+^2)]^'^'^') 

■r l(l),{2) 



-I [Sll/Ul/i2 - S12 



S22lm 



Ml M2 



|(1),{2) 
(1),(2) 



(81) 



The Fourier transforms of the symmetric and skew-symmetric weight functions are defined re- 
spectively as the jump and the average of the singular displacements across the plane containing 
the crack: 

[U](e) = U+(0 - U-(0 = ^{i sign(0 Im(Y(i) - Y^^)) - Re(Y(^) + Y^^'^)]± 
(U)(e) = \ (u+(0 + U-(e)) = ^{i sign(e) Im(Y(i) + Y(2)) - Re(Y(i) - Y(2))}s 

(82) 

We have finally recovered the expressions (j34p and (|35p , previously derived from the direct solu- 
tion of the Riemann-Hilbert problem for the interfacial crack by means of the Stroh formalism, 
consequently, we can say that the two alternative formulations are perfectly equivalent, and that 
our result is proved by this further test. 
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